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Their intrinsic properties render single quantum systems as ideal tools for quantum 
enhanced sensing and microscopy. As an additional benefit, their size is typically 
on an atomic scale which enables sensing with very high spatial resolution. Here, 
we report on utilizing a single nitrogen vacancy center in nanodiamond for perform- 
ing three-dimensional scanning-probe fluorescence lifetime imaging microscopy. By 
measuring changes of the single emitter's lifetime information on the local density of 
optical states is acquired at the nanoscale. This technique to gather information on 
the local density of optical states is important for the understanding of fundamen- 
tal quantum optical processes as well as for the engineering of novel photonic and 
plasmonic devices. 
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Quantum-enhanced sensing has become one of the major fields in quantum technology. 
The quantum properties, on the one hand, allow for measuring with a precision overcoming 
the classical limitP. On the other hand, single quantum systems represent ideal sensing 
probes. Their unique properties do not necessarily utilize intrinsic quantum effects, but 
rather the fact that quantum systems are typically very small. Therefore, they can probe 
the rather large local fields of single atoms or a local environment consisting of only a few 
molecules. Optical quantum probes provide the additional advantage of reliable initialization 
as well as efficient and easy read-out. In this paper we report for the first time on utilizing a 
single quantum emitter (QE) for performing three-dimensional scanning-probe fluorescence 
lifetime imaging microscopy (FLIM) introducing the technique of QE-FLIM. 

First optical scanning quantum probes have appeared more than 10 years agcP^ in flu- 
orescence microscopy. By now, recent advances with respect to the stability of solid-state 
emitters have paved the way for several novel approaches where quantum emitter scanning 
probes are used to detect electrkPor magnetic fields 6 . It is also possible to detect vacuum 
fields via modifications of the spontaneous lifetime of an emitter, since it is known from 
quantum optics that the spontaneous emission is not an intrinsic property of an emitter, 
but is determined by the local density of states (LDOS). More generally speaking, any light- 
matter coupling, as described on the fundamental level by a single dipole coupled to modes 
of the electromagnetic field, can be modified by changing the LD09^. By carefully designing 
the LDOS, it is possible to significantly enhance the functionality of devices in photonics 
and plasmonics. Examples are spontaneous emission control for fast optical modulator^, 
for energy-efficient lasing 9 , or for improving light trapping in solar cells 10 . The design of 
the LDOS is particularly crucial for fundamental few-photon devices in quantum optical 
engineering, e.g. for efficient and fast single photon sources needed in optical quantum 
computing^. Photonic structures such as microcavities 12 , optical antenna^, and photonic 
metamaterial^l allow for a design of the LDOS in all three spatial dimensions. Therefore, 
techniques to obtain precise information about the LDOS on the nanoscale are crucial. There 
exist several approaches to gain this information, e.g. coating of the structures of interest 
with fluorescent dye^, mapping with scanning near field microscope ^ 16 1 17 1 , nanopositioning 
of defect center^^ or colloidal quantum dot^Sl, or employing a scanning electron micro- 
scope^. Fluorescent probes indicate the LDOS via the observed lifetime changes, with most 
previous probes utilizing large ensembles of emitters like molecules in nanobead^ZI. However, 
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because of averaging over an ensemble with different spatial positions and electromagnetic 
environments, the excitation decay curve is multi-exponential. This makes it difficult to 
quantify modifications of the decay dynamics when scanning the probe. Additionally, the 
ultimate spatial resolution is still given by the diameter of the doped beads. 

A fundamental FLIM probe would consist of a single atom, and indeed two groups have 
shown the capability of single trapped ions to map the LDOS 4 23 . However, single ions in an 
ultra-high vacuum environment do not meet the requirement for a robust scanning-probe 
where a point-like fluorescent dipole is located at a scanning-tip, which can be actively 
stabilized and scanned across an arbitrary substrate. Defect centers in nanodiamondd^^EH 
provide optimum properties for this purpose, since they are optically stable even at room 
temperature. Their optical stability even allows for trapping with optical tweezer^SEH Here, 
we use a single-photon emitting NV center as a scanning probe for fluorescence lifetime 
imaging. The single emitter nature of the NV center leads to an increase in obtainable 
resolution, only limited by the size of a single NV center and the mechanical stability of the 
microscope. With only one emitter participating, there is no temporal broadening due to 
lifetime variations between different emitters in an ensemble. Furthermore we introduce a 
novel mapping technique to gather three-dimensional (3D) lifetime information while keeping 
the orientation of the probe fixed. The power of this approach is demonstrated by mapping 
the local density of states in the vicinity of silver nanowires^. 

RESULTS 

The experimental setup is a confocal microscope combined with an atomic force mi- 
croscope for simultaneous measurements (see Methods and Figure [5] for details). A very 
important task in quantitative QE-FLIM is the characterization of the probe. It has to be a 
single emitter with known properties such as orientation, lifetime, and quantum efficiency. 
To prove the single emitter character, measurements of the autocorrelation function g^(r) 
are performed before and after the nanodiamonds (typical sizes ca. 30 nm) are attached 
to the tip (see Figure [l] (a,b)). A g^(0) < 0.5 indicates that the main photon contri- 
bution is from a single emitter^. Only in this way averaging over several emitters with 
slightly different lifetimes and slightly different positions can be excluded and smearing out 
the measurement signal is avoided. A change in the bunching behavior when comparing 
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FIG. 1. Probe characterization, (a) Autocorrelation function g( 2 \r) of an NV center in 
a nanodiamond used in the scans lying on a glass coverslip. The red line is a fit to the data 
according to Jelezko et alpS. It yields g^ 2 \0) = 0.25. (b) Autocorrelation function of the NV 
center after being glued to the silicon tip. The red line again is a fit. It yields g( 2 )(0) = 0.31. Note 
the change in the NV center's bunching behavior as its environment changes, (c) Lifetime curve 
of the nanodiamond at the tip. The black curve is for the nanodiamond prior to glueing it to the 
cantilever. A fit yields a lifetime of 43.6 ns. The blue and red curve are for the diamond on the 
tip at the glass interface and at a silver nanowire, respectively. Lifetimes are 19.8 ns and 9.1 ns. 
Only photons after the dotted line are used for lifetime fits. (d,e) Lifetime versus height of the 
probe over a glass surface measured for two different nanodiamonds. The red curves are calculated 
according to the theory given in references! f or a dipole over a surface. 

Figures [j](a) and (b) is attributed to a changed dielectric as well as a changed chemical 
environment, which can influence the charge of the NV center^5Eil32j < The lifetime can be 
extracted from Figure [jjc). It shows a mono-exponential decay curve of the NV center's 
excited state for different situations: black for the nanodiamons on a coverslip (r = 46.3 ns), 
blue for the nanodiamond glued to the tip far from the surface (r = 19.8 ns), and red with 
the tip close to a silver nanowire (r = 9.1 ns). After the nanodiamond is glued to the tip, 
the lifetime only varies when the environment, and therefore the LDOS, changes, e.g. when 
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approaching a silver nanowire (red curve in Figure [jjc)). Figure [j](d,e) show detailed data of 
nanodiamond probes approaching a glass surface. When the emitter approaches the surface 
an oscillatory behavior is visible as predicted by the theory for a dipole over a surfac^Sl. 
Deviations can be explained by the fact that the silicon AFM tip gives rise to an additional 
surface not covered by the theory given in Reference 30 . From the theoretical curves the NV 
center's orientation, its position within the nanodiamond, as well as the quantum efficiency 
of its optical transition can be derived in principle (see Supplementary Information). 

With a pre-characterized and scanable single quantum emitter, the LDOS can now be 
mapped on the nanoscale in all three spatial dimensions with high resolution. In order to do 
this, for each position of a scanned area the detected photons are sorted with respect to the 
actual height of the oscillating cantilever (see Methods). In this way the three dimensional 
quantum emitter fluorescent lifetime microscopy (QE-FLIM) is established. 

A first example is shown in Figure [2j where silver nanowires of approximately 100 nm 
in diameter are raster-scanned by QE-FLIM. Figure [2](a,d,g) are AFM topography scans 
showing a wire, a zoom at the wire's end, and a wire network, respectively. In Figure[2](b,e,h) 
and (c,f,i) lifetimes images are shown for the emitter being in the distant quarter of its 
oscillation and in the quarter close to the surface. The emitter's lifetime decreases close 
to the surface, as it is expected due to the higher index of refraction and therefore higher 
LDOSPSI. Close to the nanowire the additional plasmonic modes account for an even more 
reduced lifetime. Sharp features of the image of the nanowire observed when the emitter 
is close to the surface blur out when the emitter is distant. In Figure [2](g-i) a network of 
crossed wires is scanned with the QE-FLIM. Here, the cantilever's oscillation amplitude was 
set twice as large as in Figure [2](a-f). When scanning along the nearly horizontal wire the 
expected decrease of the emitter's lifetime is observed, but also an increase in lifetime is 
found in a scan parallel to it. An explanation for this behavior is that the nanodiamond's 
position is not exactly at the center of the AFM tip, so that an asymmetry arises. These 
artifacts can be avoided by taking into account topography information as demonstrated in 
the following. 

A unique feature of 3D QE-FLIM is the possibility to derive topography corrected lifetime 
images. Whereas Figure [2] displayed lateral scans for two different relative position of the 
scanning cantilever, Figure [3] shows a scan crossing a silver nanowire in a plane perpendicular 
to the sample surface. The amplitude of the cantilever's oscillation of 128 nm is divided in 
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FIG. 2. Silver nanowire images by QE-FLIM. (a) Atomic force microscope image of a silver 
nanowire of diamenter 100 nm acquired in tapping mode. (b,c) Simultaneously recorded lifetime 
images for the cantilever being in the most distant quarter of its oscillation and in the one closest 
to the sample, respectively. Scalebars in (a-c) are 1 \im. (d) Zoom into the AFM image indicated 
by the dotted white square in (a). (e,f) QE-FLIM micrographs taken with the nanaodiamond 
distant and close to the surface. It is apparent that the lifetime of the NV center is longer further 
away from the nanowire and at the same time its QE-FLIM image blurs out. Scalebars in (d-f) 
are 200 nm. In (g-i) data for a network of single nanowires are shown. At the nearly horizontal 
wire in addition to the expected decrease of the lifetime an increase is found. This is is an artifact 
stemming from topography and the nanodiamond's position on the tip, which can be corrected (see 
text). Scalebars are 500 nm. Cantilever oscillation amplitudes are 37 nm for (b,c,e,f) and 74 nm 
in (h,i). 

25 height bins. At the same time for each pixel the absolute height of the sample is acquired 
with the AFM. Therefore Figure [3] represents a lifetime image with completely known spatial 
coordinates avoiding artifacts which often appear in scanning probe images. This method of 
topography correction can be applied to any sample. It is apparent how the lifetime decreases 
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FIG. 3. X-Z QE-FLIM line scan of a silver nanowire. (a) Color coded lifetime data as a 
function of height (z) and position (x). The absolute height is corrected according to the topography 
data acquired simultaneously with the AFM. In this way, artifacts are corrected as well (see text). 
Height and position axis are scaled equally and the scalbar is 100 nm. (b) Artists view clarifying 
the geometry of silver wire and data in (a), (c) Plot of the data in (a) with the gradient of the 
measured lifetime in decay rate shown as arrows. The arrows point in the direction where the 
LDOS is highest, whereas the size of arrows indicates the absolute value of the gradient of the 
LDOS. 
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FIG. 4. QE-FLIM scans of a silver nanowires. (a,b,c,d) show slices through the 3D data 
of the nanowire network in Figure [2] at heights of 53 — 57 nm, 76 — 80 nm, 99 — 103 nm, and 
137 — 141 nm, respectively. Grey areas are inaccessible due to being below the surface while black 
areas are not reached with the set oscillation amplitude. (e,f) show X-Z lines through the data at 
the positions indicated with the dashed lines in (a), (g) shows the slices from (e,f) and the AFM 
height data in a pseudo three dimensional manner. Scalebars are 500 nm. 

when approaching the glass cover slip or the surface of the silver wire. The observed rate 
enhancement lies well in the order expected from previous experiments dealing with the 
coupling of nanodiamonds to silver nanowire^. The highest rate enhancement found is a 
factor of 6.1, when comparing the diamond on the glass cover slip to the diamond coupled to 
a silver wire and a factor of for 14.5 for a gold nanosphere (see Supplementary Information). 
As an additional example Figure [4] shows lateral slices parallel to the sample with silver 
wires for a specific height as derived from the full 3D data set. 

DISCUSSION 

From the three-dimensional topography corrected data sets it is possible to derive another 
important parameter, the gradient of the fluorescence lifetime map. Figure [3](c) shows a 
gradient QE-FLIM image of a single silver wire. It is apparent that the gradient indicates 
the position of nanometer-sized objects which modify the LDOS. This may be used to 
localize such objects even in topography- free sample, e.g., buried semiconductor quantum 
dots™ As another important feature gradient QE-FLIM provides insight in the coupling 
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mechanism of the quantum emitter probe and the sample. For example as was shown 
recently^ the lifetime-reduction due to Foerster-coupling depends in a characteristic way 
on the dimensionality of substrate. Gradient QE-FLIM images can indicate these coupling 
mechanisms directly and with high spatial resolution. 

In conclusion we have introduced a scanning probe microscope employing a single point- 
like quantum emitter as a probe for the LDOS. This technique avoids broadening of the 
measurement signal due to multiple emitters at different locations or with different lifetimes. 
By operation in tapping mode and collecting the photons tagged with the cantilever position, 
lifetime data in all three dimensions can be acquired in a single scan without topography 
artifacts. These three dimensional data sets give detailed insight to the electromagnetic 
environment at the nanoscale with sub-diffraction-limit resolution. 



9 




FIG. 5. Measurement scheme of QE-FLIM. (a) A nanodiamond containing a single NV center 
is glued to the tip of an AFM cantilever. The AFM is operated in tapping mode and the cantilever 
is positioned in the focal volume of a sample scanning confocal microscope with a pulsed laser and 
an avalanche photodiode. Cantilever oscillation, coordinates, laser timing, and photon time tags 
are recorded, (b) For each pixel the photon detection events are sorted according to the current 
height of the cantilever. The lifetime is determined for different height bins hi by histogramming 
the arrival times with respect to the laser pulse and by making a fit to the tail of the exponential 
decay. This makes it possible to gather real 3D lifetime data. 



METHODS 



Confocal microscope The optical microscope in this experiment is a homebuilt sample- 
scanning inverted confocal microscope (see Figure [5](a)). An objective with a numerical 
aperture of 1.35 (UPlanSApo60XO, Olympus) is used to excite the NV center and to collect 
the single photons. Excitation is done with a picosecond laser at a wavelength of 531 nm 
with a variable repetition rate (LDH-P-FA-530, Picoquant). A repetition rate of 10 MHz 
is used when scanning with the QE-FLIM microscope, while a rate of 80 MHz is used for 
the autocorrelation measurements. Autocorrelation measurements are performed with a 
Hanbury-Brown and Twiss setup. After spectral filtering avalanche photodiodes (SPCM- 
AQR-14, Perkin Elmer) are used to detect the single photons. 

Atomic force microscope A tip scanning atomic force microscope (NanoWizard II, 
JPK Instruments) is mounted on the confocal microscope (see Figure [5](a)). The nan- 
odiamonds (MSY 0-0.05 GAF, Microdiamant AG) are spin coated on a cleaned coverslip, 
pre-characterized using the confocal microscope and subsequently attached to the AFM 
cantilever. A technique similar to the one described in Reference^ is used to craft the di- 

10 



amonds to the tip. To make the nanodiamonds stick, the cantilevers are dipped in a drop 
of poly-L-lysine solution!^. The typical size of the nanodiamonds glued to the tip is about 
30 nm. Cantilevers are silicon tapping mode cantilevers (NSG03, NT-MDT), which proved 
to contribute only very little to background fluorescence. While acquiring image data, the 
nanodiamond which is attached to the cantilever is positioned in the confocal volume of the 
microscope and the sample is scanned with a sample scanner. 

Electronic setup and data processing Photon detection event time tags are recorded 
using a PicoHarp 300 (Picoquant) data acquisition electronic. Markers for spatial coordi- 
nates are also fed to the PicoHarp 300 and are collected within the time tag stream. To 
extract the 3D lifetime information, in addition to the current scanning position, the can- 
tilever's oscillation is also fed to the correlation electronics. With that information, the 
photon detection events at each pixel can be attributed to different height bins (see Fig- 
ure^). Time tags for the cantilevers oscillation are created from the AFM's modulation 
reference output using a Schmitt trigger and additional counting electronics, so that only 
every 4096th period is time tagged. Data in every pixel then is histogrammed with respect 
to the cantilever height and fitted with a monoexponential decay to calculate the lifetime 
(see Figure [5^b)). Events in the first 5 ns are not used in order to avoid artifacts from fast 
decaying background contributions. The line scan in Figure [3] was averaged over 50 AFM 
scans to compensate for drift. 
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SUPPLEMENTARY INFORMATION 



I. PROBE CHARACTERIZATION 

For lifetime imaging with a single quantum emitter characterization of the probe is crucial. 
It is not only necessary to prove that the nanodiamond glued to the tip contains indeed only 
a single NV center, but in addition, its quantum efficiency, position inside the host crystal, 
and orientation has also to be determined. While it is straightforward to prove single photon 
emission with a Hanbury-Brown and Twiss setup, it is more difficult to obtain information 
on quantum efficiency, position, and orientation. We utilize the model of an oscillating single 
dipole in front of an interface 1 2 to derive this information. By measuring the NV center's 
lifetime at different heights over a coverslip and comparing it to the prediction by the simple 
model it is possible to calculate its quantum efficiency, position, and orientation. 

The theoretical expression for the local density of optical states (LDOS) for a dipole in 
front of the interface between two media i is given by 

P{*AA) = P~ = ^ = cos(<^]|| +sin(0)[^] ± (1) 



k r ,o Pq Pq Pq 



with 



oo 

[^]|| = \ I - sy]e 2ik ^}ds, (2) 



OO 

P _3 
[ i% U " 2 



oo 

f S 3 „ 

/ Re{— rh 2iklZoS *} ds, 
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the Fresnel coefficients r L and r" 



^2,1+^1^,2 

the relations for the wave vectors h inside medium % 



r = — ; ; ; — , (5) 



k z , = x /kf-(kl + kl), (6) 

and the abbreviations 



s — y k x + k y /k\ , 
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FIG. 6. Decay rate of an NV center on cantilever tip versus nanocrystal-interface 
distance. The measured data are shown in black while the model following Equation [T] is shown 
in red. The values for the red curve are: k nr = 9 ms _1 , k r $ — 36.9 ms _1 , <j) — 30.7. This results 
in a quantum efficiency of 0.8. 

Sz = (1 - ST 72 - 

The LDOS p{z, A, (j)) is dependent on the dipole-interface distance z, the wavelength A, and 
the dipole's orientation 0. It is assumed that pi — [12 — 1. 

To adapt this formula for NV centers, the LDOS for two orthogonal dipoles is added 
with the assumption, that both dipoles have the same free radiative decay rate. Figure [7^i 
shows two extremal possibilities, both dipoles parallel to the surface and one parallel and 
one perpendicular. Here, we approximated the broad room temperature spectrum of the 
NV centers by a Gaussian centered at 700 nm with a standard deviation of 50 nm. 

The quantum efficiency QE of the NV center can now be calculated according to 



QE 



k r $ 



L> 1 L> ' 

ki^z^j == kyij* ~\~ k r ^ kjiy* ~\~ ky* o pi^z^ A, 



(7) 



(8) 



The values of the parameters in Equation [T] now can be adjusted to fit the data measured, 
as shown in Figure [6| The offset in the z axis from the position of the NV center inside the 
nanodiamond is set to 0, since the quality of the data used and the accuracy of the model 
is not high enough to estimate values in the order of nanometers. Deviations of data and 
theory are due to the fact, that the model describes a significantly simplified system - effects 
arising from the interfaces of the AFM tip are not covered. Nevertheless by modeling the 
experimental situation in more detail, either analytically or numerically, a better agreement 
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FIG. 7. Variation of different parameters in Equation [T} In a the orientation of the dipoles 
(0) is varied. In b two different quantum efficiencies together with the measured data are shown, 
c shows the model for four different efficiencies. 




FIG. 8. A gold sphere scanned by QE-FLIM. In a the NV center is distant to the surface 
while in b it is close. Here the maximum enhancement by a factor of 14.5 is found. Scalebars are 
100 nm. 

might be achieved. In the following we utilize the simpler model and show the influence of 
the variation of different parameters. 

Figure [7] shows the effect of a variation of three parameters of Equation [TJ 



II. QE-FLIM WITH GOLD SPHERES 

Gold nanospheres are scanned with QE-FLIM. Figure [8] shows a scan over a gold 
nanosphere of 90 nm height. The maximum rate enhancement found is a factor of 14.5 
when the nanodiamond probe is close to the sphere's surface. 
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height bins: 

a) 0-26 nm 

b) 26-51 nm 

c) 51-77 nm 

d) 77-102 nm 

e) 102-128 nm 

f) 128-153 nm 

g) 153-179 nm 

h) 179-205 nm 



El 



Lifetime /ns 



FIG. 9. A gold sphere scanned by QE-FLIM. (a-h) are topography corrected images of a gold 
sphere scanned with QE-FLIM. The same data set as in Figure [8] is used. Scalebars are 100 nm. 
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